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Abstract 



When the phase structure of the black brane in grand canonical ensemble is discussed, the bubble 
O ; phase with the same boundary data should be included in this structure. As such, the phase transitions 

among bubbles, black branes and "hot flat space" are possible, therefore giving a much enriched phase 
structure. We also argue that under some conditions, either the grand canonical ensemble itself is 
unstable or there are some unknown new phases. 
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^ ■ 1 Introduction 

Understanding the thermodynamic phase structure of black holes is helpful in learning the properties 
of space-time and quantum gravity. The phase structure of the AdS black holes can also be used 
to study the corresponding field theory within the context of AdS/CFT correspondence. The well- 
known Hawking-Page phase transition [TJ between the AdS black hole and the "hot empty AdS space" 
corresponds to the confinement-deconfinement phase transition of the large- N c J\f = 4 super Yang- 
k> | Mills at finite temperature [2]. The charged AdS black hole in the canonical ensemble also displays 

a van der Waals-Maxwell-like phase transition which can also be understood from the dual-field 
theory [31 S]. 

The Hawking-Page or van der Waals-Maxwell-like phase structure is not only present in AdS 
black holes but also exists in some other black holes. The Hawking-Page- like transition can also be 
found in the canonical ensemble between chargeless asymptotic flat black holes and "hot flat space" as 
discussed by York [5]. In the grand canonical ensemble, there can also be transitions between charged 
black holes and the hot flat space [6] . The van der Waals like phase transition was also seen in 
asymptotic flat and dS black hole [8]. Unlike the black holes in AdS space, asymptotic flat and dS 
black holes are not thermodynamically stable objects due to Hawking radiation. So in discussing the 
phase structure of these objects, the black hole must be put in thermal contact with a heat reservoir 
with fixed thermodynamic data on the boundary so that they can form a equilibrium, and then the 
phase structure could be discussed. 
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Black branes are the black hole-like solutions of low-energy effective theory of string or M-theory, 
which are also asymptotic flat. They also need to be put inside a heat reservoir in order to study 
their thermodynamic properties. In fact, they have similar phase transitions as those asymptotic flat 
black holes. In [9], the phase structure of black p-branes in D-dimension (D = d + d + 2,p = d — 1) 
space-time in canonical ensemble is discussed, where the charge q inside the cavity is fixed. In this 
ensemble, for d > 2 charged cases, there is a critical charge q c . When q < q c , there is a first-order van 
der Waals phase transition between the larger black brane solution and the smaller one at a certain 
transition temperature. However, the chargeless case and d < 2 charged cases are different from d > 2 
charged cases. There is no such van der Waals phase transition and below a minimal temperature 
T m i n , there could not exist a black brane phase. For the chargeless case the hot flat space is the stable 
phase below T m ; n , whereas in the charged cases we do not know what could exist at the present stage, 
because of the absence of the charged hot flat space. One would wonder what could fill this part of 
the phase-space. 

It is well-known that there exists a regular "bubble of nothing" or bubble [1Q|, [TT] which carries the 
same flux as the black brane. One therefore wonders whether this bubble may play some role in the 
phase structure of the black brane. Precisely based on this consideration, bubble was considered in the 
canonical ensemble in [12j and the phase structure is enriched. There can be topological transitions 
among black branes, bubbles and hot flat space under certain conditions. The bubble really fills some 
part of the phase-space though there is still some phase-space left unoccupied. 

The grand canonical ensemble of black brane system is a little different. There are no van der 
Waalsolike phase transitions as in the canonical ensemble. Since the charge is not fixed in this case, 
charged black branes can undergo Hawking-Page-like phase transition into hot flat space. In (13], 
only black branes and hot flat space are taken into account in the grand canonical ensemble. In the 
present paper, we will discuss the role played by the bubble in the phase structure. The bubble can 
be generated from the black brane solution by a double wick rotation in Minkowski space-time, and 
in Euclidean space-time the effect is to interchange one world volume spatial coordinate and the time 
coordinate x O t in the metric, while leaving the form field and dilaton unchanged. So the bubble 
geometry is regular and has no horizon, hence no entropy. One expects that the period of the space 
direction L plays a similar role in the bubble case as (3, the inverse temperature T, plays in the black 
brane system. As a result, the phase structure depends not only on the temperature and potential 

but also on L as in the canonical ensemble in [12]. We first follow closely the papers [13] and [12] . 
However, we will find out that in some conditions there could be phase transition processes where the 
Gibbs energy itself is not continuous and the system tends to the boundary of the cavity. Since our 
"zero-loop" approximation breaks down when the horizon tends to the boundary in these conditions, 
this may not be a right behavior of the system. Nevertheless, if we take these cases seriously, we will 
argue that there may be a new phase near the boundary developed either by quantum effects or other 
nonperturbative effects, or the grand canonical ensemble is not stable under certain circumstances. 

In [12], extremal branes are also considered in the canonical ensemble and fill some part of the 
phase-space of d = 2, 1 cases. Extremal branes are very different from the nonextremal ones, as 
discussed in [14^ [To] I16pl . Their thermodynamic properties may not be obtained simply from the 
r + = r_ limit of the nonextremal ones. For example, they can form equilibria with the environments 
at arbitrary temperature. This is because their topology near the horizon is different from the nonex- 
tremal ones. If we only consider the Euclidean time direction and the radius direction, the topology 

4 |14[ H5| [T6] argue that even when the area of the horizon for the extremal black hole is not zero, the entropy still 
vanishes. This seems to contradict with the string theory calculations. Although there are some efferts trying to resolve this 
discrepancy, for example [T7] and the references therein, the reason for this discrepancy is still not fully understood. However 
we do not need this result in our paper since the area of the extremal black brane discussed in this paper is zero and hence 
it has no entropy. 
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of the nonextremal brane is R 2 near the horizon while the one for the extremal brane is S 1 x R. So 
there is no need to avoid the conical singularity for the extremal brane and it can have arbitrary 
period in the Euclidean time direction hence arbitrary temperature. There is also other evidence that 
the extremal black branes can not be seen as the continuous extremal limit of the nonextremal ones, 
e.g. [17] . The extremal black branes also have zero horizon area hence no entropy |18j . Because of 
their BPS nature and zero entropy property, they may also split into smaller extremal branes with 
fewer charges without changing the entropy. However, if we do not care about this, we could still 
include the nondilatonic extremal brane in our discussion. We will see that only when the potential 
conjugate to the charge at the boundary is fixed at a certain value can it exist. In this case, it has 
zero Gibbs energy which means that it can coexist with the hot flat space. For dilatonic extremal 
branes, since the horizon is already singular, our zero-loop approximation is not appropriate near the 
horizon. Therefore we would not consider them in our discussion. 

The organization of this paper is as follows: In Sec. [21 we review the reduced action or the Gibbs 
free energy for black brane and obtain the one for the bubble. In Sec. El we discuss the phase structure 
of this grand canonical ensemble including black branes, bubbles and the hot flat space. The final 
section contains some discussion on the extremal cases and the problems we meet. 



2 The action 

Let us recall the black p- brane solution in Euclidean signature in space-time dimension D = d + d + 2, 
(d = p+l) 

d d d P a 2 a 2 

ds 2 = A + A_^ dt 2 + A^{dx 1 f + A 1 ^ J^(dx i ) 2 + Al 1 Ap~ 1 dp 2 + p 2 ApdQ 



i=2 



A, ,1 _,-p Q 0o/ 2 



\ d/2 / \ d/2 



r+J VP 



-.2 



dt A dx 1 A ... A dx % 



F[ P +2] = dA [p+1] = -ie a ^ o/2 d ' dp Adt A dx 1 A ... A dx p , 

X<t>-<M 



where a is the dilaton coupling and in supergravity theories with maximal supersymmetry, a 2 = 



j^2- A± is defined as A± = 1 — (jf-^J f° r r + > r -> with r± being the two parameters 
characterizing the solution and related to the charge and the mass of the black brane. <pQ is the 
asymptotic value of the dilaton at infinity which is related to the string coupling g s = e^°. We have 
extracted one spatial direction x 1 from the sum of the spatial directions. See [191 120] for details of the 

solution. From the metric, we see that the physical radius of the d+l sphere should be p = At d p 

2 

and we also define f± = A^_ d r±. The charge is calculated to be Qd = 4=^ e~ a ^ / 2 {r + r^) d / 2 . 

As mentioned in the introduction, the metric of the bubble is obtained from the black brane metric 
just by interchanging the time coordinate t and one of the world volume spatial coordinate x 1 

d d d P a 2 a 2 

ds 2 = A~ 2 dt 2 + A+Aj^tdx 1 ) 2 + A^ ^{dx 1 ) 2 + A^A^dp 2 + p 2 A™dn 2 ~ +r (2) 

i=2 

The form field -Afp+i] and dilaton <f> are kept unchanged. The space-time is restricted to the p > r + 
region. The singularity of the form field at p = is excluded from the geometry and the charge is 
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the result of the flux around the noncontractible S d+1 . The period of the x 1 direction is chosen to 
avoid the conical singularity at p = r + and so the solution is regular. This solution can be obtained 
independently without referring to the black brane solution by directly solving the classical equation 
of motion. It can have the same boundary condition as the black brane when we put it in to a cavity. 
Since for an ensemble we only fix the boundary data, any classical solution that satisfies the boundary 
conditions can form an equilibrium with the cavity; therefore one can not exclude the bubble state 
from the black brane phase structure. 

While there is an event horizon at p = r + for the black brane, there is no horizon for bubbles. For 
black p-branes, the inverse of the local temperature is fixed to be 

Pip) = A y 2 Ai^/3* = Ay 2 Ai i/a ^± 1 1 - ^ i 

d \ r+J 

whereas the local radius of the x 1 direction is arbitrary. However, for bubbles, the inverse of the local 
temperature is arbitrary but the local period of x 1 is 

1_JL 

Up) = AfA'J^L* = A f AI 1/J ^± ( 1 - 2 (4) 

d \ f\j 

to avoid the conical singularity. /3* and L* in the above equations are the inverse temperature and 
the period seen from infinity, respectively. 

As in [9j [13], we put the black brane or bubble inside a cavity with a fixed radius p = ps- To 
establish a grand canonical ensemble, we then fix all the local quantities at the wall of the cavity: the 
inverse temperature /3, local period L in x 1 direction, local volume V p -\ of the x % {% = 2, . . . ,p) direc- 
tions, dilaton cps, and the potential $ conjugate to the charge at the boundary ps- The charge/flux 

can be now expressed using the boundary data as Qd = e~ a ^ B ^ 2 (f + f_)^/ 2 with f± evaluated 

at the boundary. We can also define the potential <3? in the local inertial frame using the form field 
= —iV^K&dt A dx 1 . . . dx p where (i, x l , . . . , x p ) are the coordinates in the local inertial frame. 
So $ is the conjugate potential for Q and one can easily obtain 



d . .1 



(5) 

P=Pb 



With this setup for the grand canonical ensemble, the classical Euclidean action for the black brane 
is obtained in [S] 



rbrane _ P LV P~^d+l -d 



/a \l/2 

(J+2)f^±J +J(A+A_) 1 / 2 -2(J+l) 




A - d l 1 "" 

where we have expressed the volume V p = LV p -\. According to the zero- loop approximation of the 
path integral, the Gibbs free energy G can be obtained from the classical action by G = Ie/P |21j . 
Since G = E — TS — V p Qd<&, we can identify the energy E as the first term in ([6]) divided by /3, which 
is consistent with the ADM mass as ps — > oo |22j . and the entropy S from the second term which is 
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the same as the one obtained in [T8]. Since the bubble metric is just obtained by interchanging t and 
x 1 coordinate, we expect to obtain the bubble action by interchanging f3 and L in the brane action: 



rbubble_ PLVp-ltofol J 

1e ~ 2^ Pb 



/A \ 1/2 

(J+2)(^±j +J(A+A_) 1 / 2 -2(J+l) 



i i i 



2k2 . + a_* -m^Q3 (7) 



As in the black brane case, one could check that by requiring the local minimum of the action with 
respect to charge Qd and r+, we recover the equation of state L = L(ps) and 3> = &(p~b) for 
equilibrium. This justifies the validity of this bubble action. Because of the absence of a horizon, the 
entropy of the bubble is zero and hence the sum of the first and the second term divided by (3 can 
be identified as the energy of the bubble, which can be checked to be consistent with the definition of 
ADM mass as ps —> oo. 

For future convenience, we define the dimensionless variables 

^V, b = -^, R=j 1 -, q =(^Y ,v = V2 K e~^/H. (8) 
p B J 4-irpB 4:irp B \PB ) 



where 



We use the barred variables b, (p, R to denote the quantities fixed on the boundary. With these 
dimensionless variables, one can rewrite the conditions for equilibrium as 

b = b(x,q), b&q)^-,- X^-> — (10) 



i_ 4) 2d [i-^v 



for black branes, while for bubbles 



1 x 1 ^(l - x) 1 / 2 
R = R(x,q), R(x,q) = 1 r , (11) 



d 



and 



ip = cp(x,q), (p(x,q) = -[- (12) 

for both branes and bubbles. We first constrain the range of x, q to < q < x < 1, and hence 
< (p < 1. In fact, this range can be extend to x = q < 1, <p = 1 continuously, which will be 
considered in the discussion section. The x = q = 1 is not well-defined and in this limit tp could have 
arbitary value. We will argue that this case can not be physically reached exactly and near this point 
our zero-loop approximation is not applicable. 
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We can also define the reduced action for black branes 

2k 2 I e 



(4vr)2^+ 2 y p _ 1 ^- +1 



(d + 2 




+ f J(l-^)V2 h_«_ 



2\ 1/2 



-bR 



and by exchanging b and R we obtain the reduced action for bubbles 

1/2 



2(d + 1) + dq(p 



(13) 



I E ubhlc (b,R,^q,x) 



bR 



(d + 2) 



te 1+1 /<* 



2 \ 1/2+1/d 



2(cZ+ 1) + dg<£ 

(14) 



The grand potential or the Gibbs free energy is proportional to the reduced action divided by b. 
Therefore, for a fixed inverse temperature b, they differ only by a positive constant factor. So, to 
compare the Gibbs free energies of the bubble and the black brane at a fixed temperature, we can just 
use the reduced actions instead. Notice that even though we use the same x in the bubble and brane 
action, this does not mean that they are equal. They are independent variables since the physical 
meanings of r + for the black brane and the bubble are different. We also do not differentiate these 
two variables in the following discussion and this could be easily understood from the context. 



3 Phase structure 

If (|10p . ([lip and (|12p have solutions for fixed b, R, (p, there could be locally stable black branes or 
bubbles in the grand canonical ensemble. To achieve this, first we can solve Eq. ([12D to obtain 



q 2 ip 2 



x 2 1 — (1 — (p 2 )x 
and by substituting (fT5|) into b(x,q), R(x,q) in (fTUI) and ([TT1) . we find 



(15) 



b(p(x) = b(x,q) = U(p(x) , for black branes, (16) 

and R(p{x) = R(x,q) = U^(x) , for bubbles, (17) 



where 



i 



U<p{x) = = — • (18) 

d(l-(p 2 Y~2 

So to solve ([10p and ([lip is to solve u = U(n{x) for fixed u and (p, with u = b for black brane and u = R 
for bubble, and then from ()15p we find q for the black brane or the bubble. Note that, in general, for 
b R, the solution for the bubble and the black brane may not be the same. We use x to denote the 
solution for the black brane and y for the bubble. If x 7^ y, the bubble and the black brane do not 
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have the same charge/flux from (|15|) and one can easily find out that the one with larger x solution 
carries more charge/flux. 

The discussion on the solutions of u = U^(x) under different conditions has already been done in 
|13| for the black brane, and can be directly used in the bubble case. The results can be summarized 
as follows, 



1. For y < V 3 < 1) U(p{x) is monotonically increasing for < x < 1. If u > U(p{\) there is no 



solution, and if < u < U^(l), there is one unstable solution. The vs x graph is shown in 



Fig. 1(a) 



2. For ip < W there is a local maximum for U^(x) with 



at 



Ur, 



< 



2 + d, 
1 ( 2 



d(d + 2) (l - (p 7 



2d V2 + d 



(19) 



(2 + d)(l-¥> 2 ) 



(20) 



(b) 



(a) When — - < (p < \ — the Up vs x graph is shown in Fig. 1(b) 

i. < u < U^(l), there is one unstable solution. 

ii. U(p{\) < u < u max , there are two solutions: the smaller is unstable, and the larger is 
locally stable with Ie > 0. 

(4(<f+i))5 



For ip < j^j, we define x g 



(d+tni-^) where u ° 



Ua (Xg 



0. The Ucp vs x graph for this case is shown in Fig. l(c 
i. < u < U&(1 
For u 



and I(x r , 



n. 



m. 



(pyj-j, there is one unstable solution with Ie > 0. 
g < u < u max , there are two solutions: the smaller is unstable, and the larger is 
locally stable with Ie > 0. 

For Utp(l) < u < u g , there are two solutions: the smaller is unstable, and the larger is 
locally stable with Ie < 0. 



By changing u into b or R, we can obtain the property of the black brane or the bubble respectively. 
Whether the final state of the system is a black brane, a bubble or the hot flat space depends on 
which one has the smallest reduced action. There are four different cases. The first one, if both b or 



R are chosen in cases [TJ 2(a)i, 2(a)ii, 2(b)i 2(b)ii, both black brane and bubble are either unstable 
or locally stable with Ie > 0. Since the hot flat space is an equilibrium state with Ie = 0, and we 
suppose there is no other unknown equilibrium with Ie < in these cases, the system will tend to 
the hot flat space, in line with the results of [13] . However, there is a problem in this claim, because 
the hot flat space may not be the global minimum of the Gibbs energy. We will come to this point in 
the discussion section. The second one, if b is chosen in cases 2(b)i or 2(b)ii but R is chosen in case 



2(b)iii, the locally stable bubble will have Ie < and the black brane is either unstable or has Ie > 0. 
Thus, the final state will be the bubble. There could also be a problem in this claim which will be 
discussed later. The third case is interchanging b and R, bubble and black brane in the previous case. 
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The final complicated case is when both b and R are chosen in case 2(b)iii where both bubble and 
black brane have negative reduced actions. In this case we have to compare these two negative ones. 

For this purpose, we can express the reduced on-shell action for the black brane and the bubble 
in the same form: for the black brane 



while for the bubble 



fbrano 



fbubble 



-bRF^x) 



where 



d 



2(d+l) 



(21) 



(22) 



(23) 



with < ip < d/(d + 2) and < z < 1, x and y being the solutions for ()10p and (jlip for the black 
brane and the bubble respectively. Fa{z) has a minimum at z m \ 



in (|20p . and is increasing in the interval z m ; n < z < 1 while decreasing in < z < z„ 
Fig. [5J As a result, of the bubble and the black brane, the larger one will have smaller action. Since 



- , which is equal to x max 
n as indicated in 



the monotonic property of b vs x or R vs x is depicted in the same diagram in Fig. 1(c) for fixed (p, 
which is decreasing for x > x max = z m i n , the black brane will be larger, i.e. x > y, when R > b, and 
inversely, the bubble will be larger when b > R. Therefore, when R > b the black brane will have 
smaller reduced action and when R <b, the reduced action for the bubble is smaller. 

Finally, the possible phase structure of this grand canonical ensemble can be summarized as follows, 

1. When < tp, the final state is the hot flat space. 

2. When < (p < 

(a) When b G (0, 17^(1)) U (u g ,+oo) and R G (0, £7^(1)) U (u g ,+oo), the final state is the hot 



flat space, 
(b) When b G (1^(1) 



and R G (0, L^(l)) U (u g ,+oo) , the final state is the black brane. 



Reciprocally, when b G (0,L^(1)) U (u g ,+oo) and R G (U^(l) 
bubble. 



the final state is the 
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(c) When U^(l) < b < R < u g the final state is the black brane. Inversely, when U^(l) < R < 
b < u g , the final state is the bubble. 

(d) There could be some coexisting states: 

• Three components coexisting point: b = R = u g . 

• The bubble and black brane coexisting phase: U^(l) <b = R<u g . 

• The black brane and hot flat space coexisting phase: b = u g , for R > u g or R £ 
(0, t/^(l)). And the bubble and hot flat space coexisting phase: R = u g , for b > u g or 

be (0,^(1)). 
4 Discussion 

In the previous section we have discussed the possible phase structure of the grand canonical ensemble 
for the bubble, black brane, hot flat space system. Compared with the system with only black brane 
and hot flat space, the phase structure depends not only on the potential $ (or the dimensionless (p), 
temperature T (or inverse temperature j3, dimensionless b) but also on the size of one world volume 
spatial direction R. However there could be some region of the phase-space where our zero-loop 
approximation fails and we must exclude this region from our discussion. To see this, let us recount 
the phase structure here in another way. If we fix d/(d + 2) < (p, there is nothing interesting but the 
hot flat space. If we fix < (p < d/(d + 2), there could be different phase transition processes: 

1. If we fix the temperature such that b £ (0, t/^(l)) U (u g ,+oo) and increase the size of the x 1 
direction from small to large, at first the system will be the chargeless hot flat space. Then 
just above R = U<p(l) there is a phase transition to bubble with f + close to the boundary of 
the cavity and the charge close to the extreme value. Put another way, the bubble suddenly 
appears from the boundary of the cavity. This is an instantly charging process with the Gibbs 
free energy jumping discontinuously from zero to a negative one, which characterizes a "zeroth- 
order" phase transition. Then as R increases, r + and the charge are decreasing gradually while 
the Gibbs energy is increasing. After R increases to R = u g , where f + and Qd decrease to 
the corresponding values at x g and the Gibbs energy increases to zero, the system begins to 
discharge and transforms to the hot flat space through a coexisting phase of the bubble and the 
"hot flat space" . This is a first order transition with continuous Gibbs energy. Then the system 
will stay in the "hot flat space". 

2. If we fix the temperature such that b £ (U^(l),u g ) and do the same thing, at first the system is 
a black brane, of which the horizon size and the charge do not change as R increases. There will 
be a phase transition to bubble at just above R = U^(l) where the charge suddenly increases to 
near extremal one and the r + also suddenly rises to near boundary. In other words, the same as 
in the previous case, the near extremal bubble appears from the boundary. Similar to the process 
in the previous case, the Gibbs energy decreases discontinuously which is a "zeroth-order" phase 
transition. Then, as we increase R, the bubble becomes smaller and discharges gradually. When 
R is raised to R = b, where the size of the bubble has shrunk to the same size as the original 
black brane, the bubble begins to transform back to the black brane through their coexisting 
phase. This is also a first order transition. After that, the system remains in the same black 
brane phase, i.e. the charge and the horizon size no longer change. 

We can also fix R and raise b as in above two processes, and the consequence is that the black brane 
and the bubble will exchange their roles. 

We have seen from above two cases that there are "zeroth-order" phase transition processes where 
the Gibbs free energies are not continuous. In fact, in these cases, for R < U^(l), the Gibbs free energy 
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around x = 1 for the bubble, although not being at the stationary points, have already dropped below 



where only the "hot flat space" exists, as the temperature is raised to a certain value, the reduced 
action near x = 1, still not being at the stationary points, will also drop below zero. It seems that, 
if the Gibbs energy here is correct for this thermodynamic system in these situations, the "hot flat 
space" is not a global minimum of the Gibbs free energy, and the system tends to x = 1 bubbles or 
black branes which are not the stationary points of the Gibbs free energy at first sight. Nevertheless, 
from (|15p . there seems to be an extremal solution for x = q = 1 with arbitrary <p < 1. Since extremal 
branes can form equilibria with the environments having arbitrary temperature and R [141 [Pol I16j . 
one may think that this extremal brane (from the extremal limit of (pQ), there is no difference between 
the bubble and the brane) will be the final steady state in all these cases. However, this is a spurious 
solution of (|15j) . To see this, we define r = {p/pB) d an d look at p at arbitrary p which is just ([5]) with 
Pb changed to p, and using ([15]) . we find: 



From this equation, when r + = pB which is the extremal case with x = 1, the only solution has 
p(r — > 1,1) = 1 and at the equilibrium, (p = 1, which means that the x = q = 1 extremal one can 
only have (p = 1 boundary condition in equilibrium. This contradicts with the previous arbitary (p 
solution. In consequence, the solution with x = q = 1 and arbitrary (p could not be physically realized 
exactly. If we first set r = 1 exactly (which is p = pb) and x < 1, i.e. we set the boundary condition 
first and then choose f + < p~B, we find the boundary condition p(l, x) = (p which could be arbitrary. 
This implies that the x = q = 1 solution with arbitrary (p can only be seen as a limiting case which 
can not be reached exactly. This argument is similar to the one used in [23]. So we should exclude 
this phase from the phase structure. Let us look at what happens near this extremal case but not 
reaching it exactly. For black branes, as x and q are approaching 1, the curvature singularity at r_ is 
coming closer to the horizon and also closer to the wall of the cavity. Since the quantum effect will be 
essential near the singularity, our zero- loop approximation is not applicable in these situations. The 
Gibbs free energy near x = 1 may be modified by quantum corrections. Similarly, for bubbles, when 
x, q — > 1 quantum effects must be important near p = r + and the wall. So, these zeroth-order phase 
transitions may not be the correct behavior of the system and are just indications of the failure of 
our method in describing the system under these circumstances. One possibility of the system near 
x = 1 is that there is a new stationary point near x = 1 developed by quantum effects or by some 
other unknown topological solutions, which may lift the Gibbs free energy near x = 1. If the new 
stationary point is a global minimum, there will be a new phase. The other possibility is that there 
are no new stationary points and the system itself is not well-defined or is unstable, which means we 
can not set up a steady grand canonical ensemble under such circumstances. 

It is worth mentioning that some solutions with both bubbles and black holes present were found 
in five dimensinal pure gravity |24[ [25] and can be uplifted to string theory by embedding in 11 
dimensions and a series of dualities |26[ [27] . There could be configurations of branes connected by 
KK bubbles. It will be interesting to investigate whether the solutions obtained this way could play 
a role in the phase-space of the black brane and may also be related to the issue raised above. We 
will leave this possibility as a future research direction. 

After excluding the x = q = 1 extremal case, we can also study the extremal nondilatonic brane 
with x = q < 1 whose geometry at the horizon is not singular, for example, M2 brane and Mb branes 
in M-theory, and D3 brane in string theory. From (I13D . we see that it has zero entropy automatically 
because of its zero horizon area which is already known in [18] , From (|15p we obtain (p = 1 and hence 
the Gibbs energy is zero for any b from (j!3p . Thus, it has the same Gibbs free energy as the hot flat 



zero or below the one for the original black brane. Similarly, in the cases when we fix (p 6 (d/ (d+2), 1) 




ip(r — x) 1 / 2 



(24) 
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space. Therefore, only at (p = 1 could the extremal brane exist and it can also coexist with the hot 
flat space. We should point out that this is a zero-loop result. If the higher-order fluctuations are 
taken into account, the Gibbs energy of the extremal brane and the hot flat space may be different 
and these two states may not coexist. Since in our paper, we are considering only the leading-order 
approximation, we will list the coexisting state as a possible phase. As is mentioned in the introduction 
section that the extremal brane can not be seen as a continuous extremal limit of the nonextremal 
one, it must appear as a result of some noncontinuous processes such as pair production or quantum 
tunneling. But the dynamics behind these processes is not clear to us and is beyond the scope of this 
paper. As for the extremal dilatonic brane, since the space-time geometry at the horizon is singular, 
our method does not apply here and hence we exclude them from the phase structure. 

So, if we take into account the discussion in this section, the phase structure listed at the end of 
the previous section would be modified. First, another extremal brane and hot flat space coexisting 
phase with <p = 1 should be added to the phase structure. Second, in some regions of the phase-space, 
the final states of the system is not clear to us due to the limitation of the method and we should 
exclude these regions from our discussion. In particular, when (p > and either b or R is less than 
some value about which the zero-loop approximation is not applicable, we are not sure what happens. 
We could estimate this value to be (1 - <p 2 ) 1 / 2+1 / d / \2{d + 1)(1 - (p)) by requiring the reduced action 
to be zero at x = 1. For < p < d/(d + 2), we are still not sure about the phases when either b or R 
is in (0, U(p(l)). We can not say too much about the final phases in these parts of phase-space at the 
present stage. This problem could be left for future study. 
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